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Abstract This paper is devoted to the analysis of some fundamental problems of linear elasticity in ID 
continua with self-similar interparticle interactions. We introduce a self-similar continuous field approach 
where the self-similarity is reflected by equations of motion which are spatially non-local convolutions with 
power- function kernels (fractional integrals). We obtain closed- form expressions for the static displacement 
Green's function due to a unit (5-force. In the dynamic framework we derive the solution of the Cauchy 
problem and the retarded Green's function. We deduce the distribution of a self-similar variant of diffusion 
problem with Levi-stable distributions as solutions with infinite mean fluctuations describing the statistics 
Levi-flights. We deduce a hierarchy of solutions for the self-similar Poisson's equation which we call "self- 
similar potentials" . These non-local singular potentials are in a sense self-similar analogues to the ID-Dirac's 
5-function. The approach can be the starting point to tackle a variety of scale invariant interdisciplinary 
problems. 
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1 INTRODUCTION 



The discovery of fractal geometry in nature by Benoit Mandelbrot in the seventies of the last century had a 
great impact on science and opened us a new view on nature [1, 2, 3, 4]. In the last decade also an increasing 
interest has risen in the mechanics of materials having a microstructure which is exactly or statistically 
invariant over a wide range of scales. Such materials can be in a good approximation idealized as "fractal" 
materials. Such an idealization is to be understood in the same way as the notion of the "infinite medium". 
Despite there is no infinite medium in nature, this notion has proven to be a highly appropriate idealization. 
Some interesting results to describe the mechanical properties in fractal materials have been achieved recently 
[7, 8, 9, 10]. However a generally accepted theory of "fractal mechanics" is not yet established so far. Therefore 
simple even simplistic models capturing only a few aspects of the characteristic fractal behavior in materials 
are highly desirable. 

The point of departure in the present paper is the continuum limit of the ID medium which has been 
defined as "quasi-continuous linear chain" in [13]. Before we start with our demonstration we should like to 
specify "our" notion of "self-similarity": The notion self-similarity used in the present paper as well as in [13] 
corresponds to the notion of "self-similarity at a point" commonly used in the mathematical literature [6]. An 
object is self-similar in the strict sense if it consists of parts which are exact rescalcd copies of the entire object. 
In contrast is the notion "self-similarity at a point" (with "point"*^ we mean a space-point which is contained 
by or attached to the object, where we place for example the origin of a coordinate system): An object which 
is "self-similar at a point" contains only a single part which is a re-scaled copy of the entire object [6]. The 
notion of "self-similarity at a point" implies that the scaling invariance repeats over an infinity of scales. In 
this paper as well as in [13] we call an object just "self-similar" without always mentioning that we actually 
mean "self-similar at a point". To be more precise: When we call a function A(/i) self-similar with respect to 
h which we define by A(Nh) = N^A{h) for a prescribed scaling factor N, then strictly speaking the function 
A{h) is self-similar at the point ^ = 0. The simplest self-similar functions of this type are power-functions h^. 

The paper is organized as follows: In section 3 we deduce static Green's function of displacements due 
to a 5-unit force in closed form. In section 4 we construct the basic solution of the Cauchy problem in the 
form of two integral kernels which give the solutions due to prescribed initial displacement- and velocity-field, 
respectively, in analogy to the classical d'Alembert's solution. One of these two kernels determine also the 
retarded time-domain dynamic Green's function. 

In section 5 we analyze a self-similar diffusion problem which we formulate as an initial value problem 
for a given initial distribution by utilizing our self-similar Laplacian. The solutions of this problem turn out 
to be Levi-stable distributions with infinite mean fiuctuations. As a spin-off result we obtain a hierarchy of 
functions which we call "self-similar potentials" as they play the analogue role for our self-similar Laplacian 
as the Newtonian potentials do in the case of a "traditional" Poisson-equation. These hierarchy of potentials 
includes also the static Green's function as a solution of a ^-type source. 

2 The self-similar elastic continuum 

We consider a purely elastic medium in ID with a spatial continuous constant mass distribution of density 1 
where any spatial point x represents a material point which interacts harmonically with an ensemble of other 
material points located at x it hN^ (s € Zq) in such a way that the elastic energy-density at x is a self-similar 
function with respect to h (at point h = 0). This is expressed below by relation (3). The fractal properties of 
the dispersion relation of such a system have been considered in our recent paper [13]. As a point of departure 
we invoke some results from that paper: 
The Hamiltonian of that material system is 




(1) 



^ Where that "point" is a fixed-point of that scaUng operation. 
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where x denotes the space- and t the time coordinates. V{x, t, h) indicates the elastic energy density^ 

1 oo 

V{x,t,h) = - ^ N-^'[{u{x,t)-u{x + hN',t)f + {u{x,t)-u{x-hN',t)f\ (2) 

s=—oo 

which converges in the range < 5 < 2 and where we assume h > and N being a prescribed scahng factor. 
Without loss of generality we can restrict ourselves to A'^ > 1 {N G K). In (1), (2) u and u = stand for the 
displacement field and the velocity field, respectively. (2) has the property of being self-similar with respect 
to h at point h = 0, namely 

V{x,t,Nh) = N^V{x,t,h) (3) 

As a point of departure for the approach to be developed we evoke the continuum limit of (2) and the 

resulting equation of motion. We define the continuum limit asA^ = l-|-C(0<C<<l)sor = hA^^ becomes 
a continuous variable and we can write a self-similar function A(h) which fulfills a self-similarity condition (3) 
asymptotically as [13] 

m = N-'^f{N^h) « ^ / ^dr (4) 

s=-oo Jo T 

having the form of a power function A(/i) = const h^. Both the discrete as well as the continuous representation 
of (4) converge for sufficiently good functions (see details in [13]). From (4) follows that in that continuum 
limit we can write the elastic energy density (2) as a functional of the displacement field u{x, t) in the form 



r°° t) - u{x + T, t)f + {u{x, t) - u{x - r, t))^} 

which exists as in the discrete case (2) in the band < 5 < 2. The equation of motion (self-similar wave 
equation) has then the form [13] 



^2 

-Q^u{x, t) = A(s^h,o^ix^ (6) 

The right hand side of (6) can be conceived as the self-similar Laplacian of the medium. The self-similar 
Laplacian is necessarily a non-local and self-adjoint negative definite operator and is obtained as 

r"^ {u{x - t) -\- u{x -\- t) - 2u{x)) 

C Jo T 

and exists in the interval < 5 < 2. It might be sometimes convenient to rewrite (7) in the equivalent form 



^SACM^) = - L — ^1+^ — dr , < ^ < 2 (7) 



. , . d f°°{u{x-\-T)-u{x-T)), n r o /o^ 

^(^''^'^)"(")=c5d^yo ^ — V — 

We can express Laplacian (7), (8) in terms of Weyl-Marchaud fractional derivatives which are defined as 

[12] 

S f°° {u{x) — u{x — T 
and 

^'"W-fO^r '"'''^^"'"^'' ""' right -sided, 0<i<l (10) 

where < 5 < 1 indicates the range of existence of each fractional integral Di^u{x), respectively. T{z) 
indicates the F-function (faculty-function) [17] 



^^(^) = Tchl) r tS+^ ' left -sided, 0<6<1 (9) 



^The additional factor 1/2 in the elastic energy avoids double counting. 
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roo 

r{z + l)=:z\= / e-VMr, Re{z) > -1 (11) 

Jo 

The condition Re{z) > —1 is required that integral (11) exists. Then the Laplacian (7) has the representation^ 

A(^s,hM^) = -^-^TT^ + (-1)"'^') <x) , 0<S<1 (12) 

where the actual range of existence of (12) is < 5 < 2 when we write both terms as one single integral like 
in (7). We should like to mention that the representation (12) in terms of standard fractional integrals is not 
crucially helpful in order to solve problems. Therefore, we should rather conceive Laplacian (7), (8) itself as 
definition of a self-adjoint combination of fractional derivatives. 

The goal of this paper is to analyze some static and dynamic problems which are governed by the self- 
similar Laplacian (7), (8). The following section is devoted to deduce an explicit expression for the static 
displacement Green's function. 



3 Static Green's function 

This paragraph is devoted to deduce the static Green's function for the displacement field due to a (5-point 
force. We define this Green's function g{x) by 

where S{x) denotes Dirac's (5-function. The static displacement field u{x) due to a force density f{x) is then 
defined by 

A(5,,.,C)n(x) + /(x) = (14) 

and can be represented as a convolution 

/oo 
gix - r)/(r)dr (15) 
-oo 

with the static Green's function g(x) of (13) as convolution-kernel. By introducing the Fourier transforms 

1 

g{x) = — / 5(A;)e^*^^dfc (16) 

and 



1 



oo 



it follows from (13) that 



8{x) = —J e'^'^dk (17) 

By taking into account that e^^^ are the eigenfunctions of the A(5;j^^)-operator (7) with A(^g^f^^Qe'^'^^ = 
—uP'{k)e^^^ the dispersion relation is obtained from 

u5 POO . , 
a;2(fc) = y y^^^ + e"'''^ - 2) r-'^-^r (19) 

which was also given in [13] in the form 

^, = ^/;(l^d., 0<^<2 (20) 



^We have —1 = 6*" with (—1) * = e as we use throughout this paper for any complex number z = \z\e'^^ the 
value — 7r < ip = Arg{z) < tt for its argument ip. 
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This integral exists in the same interval < S < 2 just as the Laplacian (7). We use now the identity^ 

1 1 fOO 

\t\-^-^ = 7T-TTRe{ lim (e - irr^-H = lim ^ZTT^ei / e~'^'-'^h^ds} (21) 

which is well defined in the range < 6 < 2 of existence of the Laplacian. Replacing identity (21) in the 
integral (19) yields^ 

hS^ POO u5 

^'(^) = --- ^oTTt/ s'{6{s-\k\)-2S{s)}ds = —-^\kf, 0<S<2 (22) 

C5\ cos Jo (dl sm ^ 

which has the form [13] 

oj'^{k) = As\kf , 0<S<2 (23) 



^5 = 7^-— -^>0, 0<5<2 (24) 

C d!sm(^) 



where coefficient As is given by 

C 5!sin( 

We conclude that the positiveness of the constant > (which is determined by sin ^ > 0) of (24) due to 

the positiveness of integral (20) reflects nothing but the condition of elastic stability, i.e. a;^(A;) > for k ^ 
which is a physically necessary condition. This condition together with the condition of existence of Laplacian 
(7) constitute a physically sufficient condition which gives the 5-range of physical consistency namely the 
range < 5 < 2. We emphasize that the regularization performed in (21) would allow to obtain converging 
integrals \/6 > —1. However among these mathematically admissible (5-values only those within < (5 < 2 lead 
to physically meaningful results. This is true for all physical problems formulated with Laplacian (7). 

Explicit Green's function: 

The integral (16) assumes with (18) and (23) the form 

5(a^) = ^ Ifel-V'^-dfc (25) 

where we make use of the relation 

1 roo 

1^1"' = ZT-P^] ^ / e"'^^"'^^'"'dr} (26) 

(d — l)!cos ^ e-*-o+ Jo 

which is well defined in the range of existence < 5 < 2 of the Laplacian (7). Using this identity to evaluate 
(25) we finally arrive at 

g{x) = go \x\^-^ (27) 

with the pre-factor 

1 sin^ 

90 = J = Q<5<2 (28) 

^ 2^5(,5-l)!cosf 27r/i'5cosf ^ ' 

Expression (27) holds for < 5 < 2 which yields physical consistent behavior for the deformation -^g{x) 

having the exponent (5 — 2 < 0: \-^g{x — s- 0)| — t- oo being singular at .x = and | -^g{x — t- ±oo)| vanishing 
for \x\ oo. The consistence of expressions (27), (28) with (18) can be directly verified by performing the 
back transformation 

/oo roo 
g(x)e-'^'' dx = 2go lim Re{ / |x|^-^e-^(^-^'=) dx} (29) 
-oo e->-0+ ^0 

''Re(..) denotes the real- and Im(..) the imaginary part of a complex number (..). 
®See details on integrals of this type in appendix 8.4. 
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which can be rewritten as 



gik) = 2gor{S)cos^\k\-^ (30) 

This expression must coincide with the expression g{k) = ^^tq^ of (18) with the dispersion relation (23). This 
leads to the condition that 

i- = 25omcos^ (31) 
which is fulfilled by (28). It is possible to evaluate (25) in the range < 6 < 1 directly to arrive at 

g^x) = 1 lim Re{ r e-^^'-'-\k\-'dk} = ^i?e(e - ixf-' = ^(1 " '^j f 
ttAs e-)-0+ Jo T^Ag irAg 

Comparision of go of (32) and (28) lead to the condition that 

r(<5)r(i -s) = ^ (33) 

sinTTO 

which is known as Euler's reflection formula [17]. Consequently the correctness of expressions (27), (28) for 
the Green's function and (18) is proven. For the analysis to follow it is useful to generalize the definition of 
the T-function extending it to a < —1 in such a way that for the extended definition of a! Euler-relation (33) 
is fulfilled^ 

a! = r(« + l) = -— — , (34) 

r(— a) sm TTa 

where T{—a) is well-defined by (11) for a < —1. Definition (34) is motivated in appendix 8.4 for a < — 1. 
Equation (34) complements the usual definition (11) of the F-function by extending it to the "forbidden" 
a-range of (11) a < —1. 

We again empasize our previous observation that the Green's function (27)ff. is physically meaningful only 
in the range < 5 < 2 to obtain consistent behavior for the deformation -^g{x) behaving as |a:;|''~^, such as 
its vanishing for \x\ — > oo and its singular behavior for a; — >■ 0. 



4 Cauchy problem and dynamic Green's functions 

In this section we construct some fundamental solutions of dynamic problems: First in paragraph 4.1 we 
construct the solution of the Cauchy problem and in paragraph 4.2 the causal (retarded) time-domain Green's 
function. Further we make some brief remarks on the frequency domain representation the dynamic Green's 
function representing the solution of the self-similar Helmholtz-equation. 



4.1 Cauchy problem 

The Cauchy problem is defined as follows: Construct the displacement field ?x(x, t) solving (6) which we rewrite 
for our convenience now in the form 

Q^u{x, t) = -Cu{x, t) (35) 

with A(5 /j^) = —C denotes the self-similar Laplacian of (7). The operator C is necessarily self-adjoint^ and 

positive definite with Ce^^^ = uP'{k)e^^^ and the dispersion relation uj'^{k) > for A; / (Eq. (23) of the last 
section 3) . Since £ is a linear operator we can take advantage of this and employ in the analysis to follow the 
calculus of linear operators. 

The displacement field u{x, t) is required to fulfill the following initial condition at t = 0: 

®This definition of a! maintains the representations of regularized integrals of the form (126) for all admissible values a € R. 
^i.e. C = C+ 



6 



u{x, t = 0) = uo{x) (36) 

dt' 



with the prescribed initial displacement field uq{x) at t = 0. The velocity field ^u{x,t) is required to fulfil 



the following initial condition at t = 0: 

d 

_u{x,t = 0) = vo{x) (37) 

where vo{x) denotes the prescribed initial velocity field at t = 0. (35) with the initial conditions (36) and (37) 

represents the Cauchy problem which is defined in the domain — oo < x < oo and — oo < t < oo. In order to 
construct the unique solution of the Cauchy problem, it is convenient to write the displacement field u{x, t) 
at time t in the form 

u{x, t) = cos (£^t) uq{x) + £^5 sin (£^t) vq{x) (38) 
and as a consequence of (38) the velocity field v{x,t) = ^u{x,t) at time t 

v{x,t) = —C^ sm.{C^t) uq{x) + cos {C^t) vq{x) (39) 

Expressions (38), and (39) are defined by their power-series involving only entire powers with n = 
0, 1,2,.. G IUq of operator C = — ,^ with the identical operator = 1. 

To evaluate expression (38) it is convenient to write it in its spectral representation which is in our case 
the Fourier-representation. Introducing the Fourier transformations 

Mx) = TT / Mk)e"'''dk (40) 

1 /■oo 

(x) = — / vo{k)e''''^dk (41) 



and 



Vo 



and taking into account that 



with 



Cie''''' = uj{k)e'^'^ (42) 



f{Cn)e'^'' = f{uj{k)t)e'^'= (43) 

for any sufficiently smooth function f{ojt) = J2m=o am'^"^t"^ and where (jo{k) denotes the (positive) square-root 
of dispersion relation (23). It is now convenient to rewrite (38) in rather trivial manner 

/OO /" OO 

S{x -0^(0^^ + ^"^ sin (ch) / S{x-Ovo{0d^ (44) 
-oo J —oo 

where the operators act only on the dependence on x. By using the Fourier representation of the 5-function 
(17) we have 



Q{x, t)=C-"^ sin {Ch)5{x) = ^ T "^"(^ff)^) dfc (45) 

27r J-oo oj{k) 



and 



1 1 f°° 

-Q{x, t) = cos {C^t)6{x) = — / e'*^"^ cos {u{k)t)dk (46) 
^ 27r J-oo 



d 

at J-oo 

We can then write (44) in terms of the kernels Q{x,t) and §iQ{x,t) namely 

roo Q 



u{x, t) = j^^ -Q{x - t)uo{Od^ + Q{x - ^, t)vo{OdC (47) 
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and the velocity field is written as 



v{x,t) = ^Q{x - ^,t)uom^ + Q^Q{x - ^,t)vom^ 
Prom relations (45) and (46) we observe the initial conditions 

Q(x,t = 0) = 0, -Q{x,t = 0)=6{x), —Q(x,t = 0)=0 



(48) 



(49) 



which guarantee that (47), (48) indeed solve the Cauchy problem. The problem is solved for our medium by 
the explicit determination of (45) and (46). To this end we evaluate relations (45) and (46) which can be 
written as^ 



oc f.2n+l 



and 



with the functions qn{x) defined by 



2n 



n=0 ^ 



dt 



(50) 

(51) 



1 fOO T fOD 

qn{x) = C5{x) = — e'^''J^''{k)dk=:\\Tn- e-*^'a;2"(A;) cos dA: , n = 0, 1, 2, .. G IUq (52) 

27r J-oo e^o+ vr 7o 

The integral /q°° e~^^k°' cos kx dk exists for a > — 1 and contains the entire admissible range < (5 < 2 of the 
Laplacian (7). Hence we have 



ATI i-oo An 

qn(x)= lim / e-'^^k^'' cos kxdk = -^Re(Jns),n = 0,1,2,. .eKo,0< 6 <2 

€-*-0+ vr JO TT 



(53) 



with 



nS 



lim 



Jo 



e-'s''^ds= lim 



lim 



1 



, e->-o+ e — z|a; 

where (n(5)! = T{nS + 1) and denotes the principal value of (..). We obtain then 



■ TvSn (nS)\ 
le'— ,X7^0, n = l,2,.. GlU 



|/j.|n<5+l 



-KSix^ + iVi-), n = 



( AV: 



qn{x) = qn{\x\) = < 



^i?e(J„5) = -^sin(^)^|^, X7^0,n = l,2,..Gni, 0<5<2 



^ qo{x) = ^Re{Jo) = S{x) = , X 7^ 



(54) 



(55) 



where Ag is determined in (24). On the right-hand sides of (54), (55) we skipped a singular term which also 
occurs in the case n > at x = 0. Therefore this relation excludes x = for n > 0. We come back to this 
important point in section 6 by considering integrals of the Jq, for the entire range of their existence (a > — 1). 

The series (50), (51) together with functions (55) constitute the solution of the Cauchy problem being 
defined in the entire range < 5 < 2 of existence of Laplacian (7). It is shown in appendix 8.2 that the series 
for (50), (51) converge uniformly absolutely for all x and t. 



Where the convergence of these series will be verified in appendix 8.2. 
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4.2 Retarded time domain Green's function 

In view of the results of the last section it is only a small step to determine the retarded (or causal) dynamic 
Green's function. The dynamic Green's function is defined as the solution of the self-similar wave-equation 

((^ + e)' - A^sAo) *) = ^(^)'^(*) (56) 

with the self-similar Laplacian A^^^ ^ defined by (7). Causality means that g{x,t) must be zero before and 
nonzero only after the impact of the pulse S{x)5{t) i.e. g{x,t < 0) = and g{x,t > 0) 7^ [20, 21]. We 
introduced in (56) an infinitesimal positive damping term e 0+ which breaks the time-inversion symmetry 
of the wave operator and guarantees in this way g{x,t) being the causal solution of (56). The physical 
interpretation of g is the displacement field due to a S-type force density of the form of the right hand side of 
(56). The Green's function defined by (56) solves the dynamic problem 



.d_ 
di 



where f{x,t) is the density of external forces and u{x,t) the corresponding displacement field by the convo- 
lution 

/'oc rt 
/ g{x-^,t-T)f{^,T)dCdT (58) 
-00 J —00 

where g{x,t) indicates the Green's function defined by (56). The upper limit in the time-integration in (58) 
is due to the causality of ^(a:,*). With the Fourier-transformation 

1 f°° 

9i^,t) = — gik,t)e''^dk (59) 

and 

S{x) = — / e'^^ dk (60) 
27r J-00 

yields (56) the evolution equation for the normal amplitudes g{k, t) in the form 

l^{l + e)' + u'ik)yg{k,t)=d{t) (61) 

which is nothing but the equation for the causal Green's function of a damped harmonic oscillator with 
damping constant^ e > and eigenfrequency u}{k) which has the (causal) solution of the well-known form^'^ 

N / X smuj(k)t , , 

5 =e-^*e t ^ 62 

where 0(t) denotes the Heaviside unit-step function being defined as 0(t) = 1 for t > and 6(t) = for 
t < 0. (62) holds for any (also non-infinitesimal) damping e > where the dissipative factor e~^* can be 

skipped in the limit of an infinitesimal e > 0. The frequency u^ik) = v^l^l^ denotes the positive root of the 
dispersion relation (23) holding for < (5 < 2. In view of the relation 

C 2 smCHe'"'' = ttV^c'"^ (63) 

U!{k) 

with jO = —As^h,c we have the retarded space-time Green's function g{x, t) in the form 

g{x, t) = e""* e{t)C-^ sin {ck) 6{x) = e"'* &{t) Q{x, t) (64) 



e Q+ being positive and infinitesimal being sufficient to obtain the causal solution. 
^"See any textbook of theoretical physics or e.g. [21]. 
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where Q(x,t) has been evaluated in relations (50) fF. 



Some remarks on the dynamic Green's function in the space-frequency domain 

For the sake of completeness we give a brief representation of the frequency-space Green's function, i.e. Green's 
function of the self-similar Helmholtz equation. The Helmholtz Green's function is defined as 

/•oo 

g{x,u) = / e''''g{x,t)dt (65) 
Jo 

where g{x,t) is the retarded space-time-domain Green's function (64) and the lower integration limit is due 
to causality. (65) solves the self-similar Helmholtz-equation 

((w + ief + A^sAo) 9{x, uj) = -5{x) (66) 

where e > is a positive infinitesimal damping constant which guarantees existence of the Helmholtz Green's 
operator defined in (67) below and causality in the time domain. We can invert this Helmholtz equation with 
^(<5,h,C) — ~^ where the Helmholtz Green's function can be written as 

5(x,w) = (£- (w + k)^)^^(5(x) (67) 
with the Green's operator (resolvent operator) of the frequency-domain (£ — (a; -|- ie)^) ^ . 



5 Self-similar diffusion problem - Levi flights 

This section is devoted to the analysis of the diffusion problem governed by the self-similar Laplacian (7). 
We consider an ensemble of particles with local concentration p(x,t). We also can conceive this problem 
to be defined below as a fractional differential equation for a (probability)-distribution p{x,t). The idea 
to represent statistical distributions as solutions of fractional differential equations is actually not new, we 
refer for instance to the paper of Li et al. [14] where "stable" distributions are generated by fractional 
differential equations. However, so far the fractional differential equations generating statistical distributions 
as solutions were rather a kind of "guesswork" without having an intuitive physical "picture" about the 
generating fractional differential equation. The goal of this section is to analyze a diffusion-equation employing 
the self-similar Laplacian (7) of the form 

^^p{x,t) = -Cp{x,t) (68) 

where —£, = A(5 ;j g) indicates the self-similar Laplacian (7). We consider this problem in the domain — oo < 
X < oo and t > 0. It is important that we restrict on t > as diffussion is a time-irreversible process-^^ 
which is expressed by the odd-order of the time derivative in (68). p{x, t) denotes the density of the diffusing 
particles but can be also conceived as a probability-density distribution-function. (68) recovers the traditional 
diffusion equation by replacing the self-similar Laplacian ^^^^ by the traditional ID-Laplacian leading 
to a gaussian distribution. 

We give an intuitive physical picture of the diffusion processes described by (68) as follows: Let us assume 
we have p{x, t)dx particles at time t at location x. Any particle is allowed to jump at time t from position x 
to arrive at t + dt at any location x it r different from x (with < r < oo). The jump-rate of this process is 
measured by {p{x + T, t) + p{x — T,t) — 2p{x, t))/T~^^~^^^ which is given by the integrand of (7) where < 5 < 2. 
We observe the following: The jump-rate from x to x it r is zero if (p(x + T,t) + p{x — T,t) — 2p{x, t)) = 0, i.e. 
there are no jumps for a equal distribution p = const. Jumps to "far" locations are more seldom as jumps 
to "close" locations for equal {p{x + T,t) + p{x — T,t) — 2p{x,t)) / 0. The possibility of a particle to make 
jumps of arbitrary distances is expressed by the non-locality and continuity of Laplacian (7). Particle-jumps 
to any location are possible in (68), whereas non-local jumps are suppressed in the traditional (gaussian) case 

which is reflected by the locality of the traditional Laplacian The non-local Laplacian (7) sums up over 
^^That means the time-inverse process is physically inadmissible. 
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all possible jump-events by which particles can reach to or escape from a space-point x. For p = const the 
jumping rate is zero (since CI = 0) and, as we will show, the equal-distribution is a stationary solution and 
p = an attractor for t — t- oo for any initial distribution po{x) just as in the case of gaussian distribution. 

If we assume p{x, t) to be the probability distribution of a single particle and t)dx to be the probability 
to find a particle in the interval within [x, x + dx] at time t, then the trajectory of the particle propagating 
due to probability-distribution p[x, t) which evolves according to (68) is discontinuous and erratic whereas in 
the gaussian case the particle-motion describes a continuous trajectory. The discontinuous characteristics of 
the particle trajectory associated with (68) becomes more and more pronounced as (5 ^ which is due to 
the weaker decay of the kernel t~(^+^) for r — t- oo. Wc should like to mention that motions characterized by 
jump-probability distributions of the form of a power law as t~^^~^^^ are known in the literature as Levi-flights 
where S is called Levi-parameter. Levi flights correspond to Levi-stable probability distributions for < (5 < 2 
[3]. Typical features of Levi-distributions are their infinite variances in the range < S < 2. The present 
analysis will end up in probability-distributions of exactly this type. 

In order to define diffusion problem (68) as initial value problem, we prescribe an initial distribution po{x) 
at t = as 



p{x,t = 0) = po{x) (69) 

We assume that the particle number is conserved in time, there is neither particle generation nor an- 
nihilation. The diffusion equation (68) describes then the continuity equation accounting for the particle 
(probability) balance at space-point x and time t having the form 

^^p{x,t)+V ■j{x,t)=0 (70) 

where V = ^ denotes the usual ID gradient operator. Unlike in the Gaussian case where the traditional 
Laplacian comes into play by assuming Pick's law as a consequence of a gradient dynamics as constitutive 
law, in our case the constitutive law is instead of a (local) gradient a non-local constitutive law for the particle 
(probability) flux. An explicit form of the particle (probability) flux is obtained in view of (8) as 

,(..t)=-^r w'^+^-'>-/'---''» d.. o<*<2 (71) 

Co Jo T 

which is our self-similar constitutive law and can be conceived as the self-similar analogue to Pick's law 
jgauss{x,t) = —Vp{x,t) of the Gaussiau case. In (71) an irrelevant integration-constant has been neglected. 
We can consider (71) as definition of the self-adjoint part of the self-similar gradient operator acting on p{x, t). 
The distribution function p{x,t) is required to fulfil the normalization condition 

/oo 
Po(a;)dx = l, Vi>0 (72) 

-oo 

(72) indicates that p{x, t) has to be a normalized distribution. p{x, t) is uniquely determined by (68) and 
(69) and can be written in the form 

p{x,t) = e-^'po{x) (73) 

where we restrict us to t > 0. The question to be answered in the following is: Does (73) represent a normalized 
distribution, i.e. fulfills (72) under the condition that po{x) represents a normalized distribution? We can 
answer this question by "yes" by giving the following formal prove by integrating (73) 

/'C<j roo 
p{x, t) dx = e"^* / pq{x) dx = e"^* 1 = 1 (74) 
-oo J —oo 

where we take into account that £"1 = (5„o since Laplacian (7) applied to a constant yields zero. Let evaluate 
(73) and deduce the propagator (kernel) that represents p{x,t) in terms of a convolution 

/oo roo 
W{x - e, t)poim = e-^* / Six - OPoiO (75) 
-oo J —oo 
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The propagator W{x,t) is obtained by 



W{x,t) = e-^^S{x) (76) 

and fulfils the initial condition 

Wix,t = 0) = d{x) (77) 
The necessary and sufficient condition that W{x, t) describes a normalized probability distribution is 



f 



W{x,t)Ax = \ (78) 



which is indeed fulfilled as can be seen by putting Pq{x) = b{x) in (74). If p{x, f.) is conceived as a probability- 
distribution W{x, t) is the conditional probability density to find a particle which was at t = located at a; = 
at time t > at point x. The propagator W{x,t) can be evaluated by using jCe^''^ = (jj'^{k)e^^^ in the form 

W{x, t) = e-^*5{x) = — / e^^^e"'^ dk = - e'^'^ ^ cos {kx) dk (79) 

where oj^{k) = A\k\^ is the dispersion relation (23) which is defined in < 5 < 2. Integral (79) exists in the 
whole range of < 5 < 2 where As > ioi t > 0. We observe further because of u}^{k) > that e~'^ ^^^^ 
at t — >■ DO thus (79) tends to zero as t tends to infinity 

lim W{x, t) = (80) 



so the stationary solution poo = is always an attractor independent of the initial distribution pq{x) in the 
physical admissible interval < 5 < 2. It is important to note that for any finite <t < oo (80) maintains 
in some regions non-zero (positive) values in order to fulfil (78). The goal is now to evaluate (79). 

It can be further seen in (79) since e^"^^^ * is a monotonous decreasing function in k that the positive 
contributions due to coskx dominate the negative contributions in the last integral of (79) so that W{x,t) 
fulfills another necessary characteristics of a probability-distribution, namely 

W{x,t)>d VxGR andt>0 (81) 

So we can infer that the initial value problem defined by (68), (69) with the self-similar Laplacian (7) indeed 
describes a diffusion problem (transport problem) under the condition that (69) is a normalized distribution. 
It remains us now to evaluate (79) in more explicit form 



1 



oo 



W{x,t) = -Re^J 6*'=^-^^'= ' d/c| , 0<6<2 (82) 

where < (5 < 2 is the range of existence of As. Distributions of the form (79), (82) have already been 
proposed by Levi [2, 5] and refers to the category of a "Levi-stable" or "L-stable" distribution (for further 
details on L-stable distribution we refer e.g. to [2], p. 106ff. and the references therein). For ^ = 1 we can 
directly evaluate (82) directly and obtain 

oo 



Ws=i (x, t) = Re\- r~ e-*^(^i*-^^) dk \ (83) 



vr Jo 



which yields 



where Ai = Relation (84) represents spatially a distribution of the Cauchy type with Ait being a 
"parameter". One verifies directly that (84) fulfills the conditions (77), (78) and (80). Levi already gave static 
spatial probability-distributions of the general form (82) [3] and found that in general they are admissible for 
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< S < 2 where 6 = 1 corresponds to the Cauchy-distribution and 6 = 2 to the Gaussian distribution. For 
6 = 2 the Laplcian (7) does not exist and hence A2 is not defined in our approach. 

Mandelbrot showed that L-stablc distributions play a crucial role in the description of economic processes 
and stock courses and describe the statistics of fractal irregular trajectories. The irregularity of these tra- 
jectories decreases from 5 = to (5 2 [3] which is consistent with our above given physical picture. The 
irregularity is due to the fact that (82) produces in the admissible range of existence of the Laplacian (7) 
< 6 < 2 infinite fluctuations due to discontinuous trajectories. The Cauchy distribution {6 = 1) itself refers 
also to this category. Only the Gaussian case corresponds to a diffusion equation with a traditional Laplacian 
describing continuous trajectories (Brownian motion) with finite mean fluctuation (variance). A consequence 
is that the unjustifled use of Gaussian statistics underestimates fluctuations and risks! We refer in this context 
to [1, 2, 3, 4] and the references therein. We will reconsider briefly this important point at the end of this 
section. 

It still remains to evaluate (82) which we do by expanding e~'^ in a series and arrive at the expression 

Wix,t) = J2i-ir-qnix), 0<6<1 (85) 

77-. 

where this series converges only in the range < 6 < 1 which is shown in appendix 8.3. In (85) appear the 
functions 

qn{x) = C6{x) = — / e'^"" uj'^'\k) dk (86) 

2lT J -00 

which where already determined in the last section by (52)- (55). Due to their importance for our approach 
we devote the next section 6 to their thorough analysis. The series (85) writes 



^ ' ^ 7re->o+ n!(e-«|x|)"'5+i J ^ ^ tt 6->o+ n!(e - «|x|)"'5+i J ^ ^ 



which assumes for a; 7^ 



n ' — 1 ni 

n=l 



where the singular parts at x = vanishing for x 7^ have been omitted. (88) has the form 

W{x,t) = --l-^Im{w{^{x,t))} (89) 

when we introduce the variable 

^{x,t) = j^e'^As (90) 

with the function if (^) 

^c) = f:(-ir^r, o<6<i (91) 

n=l 

where series representations (85), (87) converge in the interval < 6 < 1 (appendix 8.3). 

Let us now discuss some of the main characteristics of the distribution W{x,t) which is defined in the 
whole range < 6 < 2: In view of (79) and (55) we observe that W{x,t) = W{—x,t) = W{\x\,t) is a 
symmetric distribution with respect to x and as a consequence all moments of odd order 

/oo 
W{x,t)x'^''+^dx = 0, n = 0,1,.. G lUo (92) 

-00 
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are vanishing especially the mean- value < x >= 0. The moments of even order are obtained from 

/oo 
W{x,t) x^^'dx^oo, n = l,..GriJ (93) 

-oo 

are diverging. We can conclude without any further calculation that (93) is divergent since any real valued 
power function (except a constant) is not admissible in (7). The existence of (7) however is achieved for 
certain complex valued power functions as we show in the next section 6. 



6 Self-similar potentials 

Finally we analyze more closely functions of the type qn{x) of (86) which were determined in (52)- (55). To this 
end let us evoke the definition of the qn{x) = £"5(x) which was defined by using the Laplacian (7), (8). We 
should emphasize that A^^ /j = —C depends on the continuous parameter 6 which was restricted in (7), (8) to 
the interval < S <2 where the upper limit S < 2 came into play since 2u{x) — u{x + r) — u{x — r) — )■ C{x)t^ 
is a quadratic function of r for small r. In this subsection wc will show that for a certain class of singular 
functions, namely functions of the form hnsix) = C^5{x)/A^ , the upper limit 5 < 2 does not exist and for which 
the condition 5 > is sufficient. We call these functions due to their importance as solutions of the Laplace- 
equation (7) as "self-similar potentials". Their singular behavior at a; = becomes especially important when 
we consider integrals or convolutions of these functions. Whereas any real-valued power-functions x^ are not 
admissible in (7). Self-similar potentials to be analyzed, however, are admissible and give the only way to 
an admissible definition of powers Let us consider the Fourier transformation of l/cl^e"'*^!^ with e > 0, 
namely^^ 

6„(x) = lim — / e^'^^-l'^l'litrdA; = lim - / e'^'' k'' cos kx dk , a > -1 G IR (94) 

which exists for a > —1. This function can be identified with the Qnix) of (86) for a = nS. We can rewrite 
(94) by introducing the integration variable s = fc(e — i|a;|) in the form 

1 1 /""^ „ ^ , 1 a! „ f ^ 



ba(x) = Re { lim , / e'^s'^ds } = —Re < lim — - } , a > -1 e W. (95) 

which holds for a > —1 and where a\ = T{a -\- 1). We note that (x -|- ie)~"~^ assumes the complex conjugate 
value when replacing x -H- —x. So its real-part depends only on |x| and we can rewrite 

a! „ f. . ^ 



baix) = ba{\x\) = -Re hm . . \ , a > -1 G IR (96) 

Expressions (95), (96) are in accordance with with the result given by Gel'fand and Shilov for the Fourier 
transform of |A;|" (page 447, equation (13) in [23]). 

We summarize the following important observations: 

(i) ba{x) = ba{—x) = bai\x\) is a symmetric function in ,x for all a > —1. 

(ii) ba=o{x) = 6{x) represents the usual Dirac-^-function, i.e. is localized at x = 0. 

(iii) For a = 2n with n = 0,1, 2, .. G Mo (95) takes the form of even derivatives of the (5-function 

b,r^{x) = i-ir^6{x) (97) 



For further details on the Fourier transform of \k\", see [23] 
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which are localized at x = 0. 



(iv) a > -1 G IR: 

For X 7^ (95), ba has the exphcit form (for x 7^ we can put directly e = in (95)) 

, , , a! /CKvr. 

^"^"^^^-^^^'^^^T^ ^^^^ 

where we directly verify that this expression is zero for a = 2n (in accordance with case (iii) for x ^ 0). For 
a > — 1 G IR ^ 0, 2,4, .. expression (98) is non-zero for x 7^ and hence in contrast to the Dirac's (5-function, 
non-local. Expression (98) includes also the case of odd integers a = Q!„ = 2n + 1 where n = 0,1, 2, .. G lUo- 

We observe that 6q(x) defined by (94) can be conceived as the Fourier transform of which exists only 
for a > —1. For exponents a < — 1 the integral (94) diverges. However, we can define a regularized Fourier 
transform in the spirit of generalized functions as given in Gel'fand and Shilov [23]) which exists also for 
a < — 1. We devote appendix 8.4 to this case. 

Let us consider now the following integral 

/■oo f , a > 

ba{x)dx = / e-l*^l'|A;|"(5(A;)dA; = < 1 , « = (99) 

[00, -1 < a < 



(v) For a > we observe that integral (99) is vanishing. On the other hand we observe that 6q,(x) 7^ 
for X 7^ and sign{ha{x)) = —sign{sm (^)) is the the same Vx 7^ 0. 

Let us raise the following question: How is it possible that the integral (99) is vanishing for a > whereas 
6q, 7^ having the same sign as —sign{sm{^)) almost everywhere, i.e. Vx 7^ 0? To give the answer, let us 
consider the integral 

Ta{a) = / ba{x)dx = -^^^ Sin (— ) , a > (100) 

Ja>o vra" 2 

where we have the property lim(j_j.o+ |Xc^(a)| — )■ 00. Let us again ask, why docs (100) not converge in the 
limiting case a — t- 0+ towards zero for a > as indicated in (99) since ba{x) = ba(—x) is an even function? 
If integral (99) indeed is vanishing, then (100) should be compensated by the integral 

Jaia) = / ba{x)dx = - ba{x)dx, a>0 (101) 

Jo 2 J-a 

since /q" ba{x)dx + 6«(x)dx = | /_!^ 6Q(x)dx = (6(x) = 6(|x|)). 

To evaluate this integral we have to use (94) or directly (96) which hold Vx including x = to arrive at 

1 POO / \ 

Ja{a) = Ihn^ - J e-^^k"" U cos kxdx) dk (102) 

which yields 

1 roo Co — 1 V r 1 
Jaia) = lim - / e-'^'k"-^ sin (Ifcla) dk = ^ '-Im \ 

since a > being non-zero we can further evaluate (103) by putting e = to arrive at 

Jaia) = / 6,(x)dx = ^ ^sin(— ) (104) 

Jo 7ra" 2 

which compensates indeed (100) so that Ja{a) +^a{a) = thus (99) indeed is fulfilled with the remarkable 
property 



(103) 
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lim i7a(a) = lim / ba{x)dx = sign{sm{—)) x oo , a > (105) 

a— >-0+ a— i>0+ Jq 2 

is diverging for a > as a~"! That means if we imagine for a moment \ba\ being a mass density, that the 
total mass concentrated in the interval < a; < a is the same as the total mass in the interval a < a; < oo. We 
conclude that 6a(x) has, in dependence on a, several oscillatoric peaks of alternate signs localized for e — >■ 0+ 
infinitely close to x = 0. The integral over these localized oscillations up to a finite limit a > compensates 
integral (100). The strength of these oscillatoric peaks are "infinitely stronger" as that one of a ^-function. 
The closed form expression of the which takes into account this singularity is therefore (95) which fulfils 
for instance the "desired" property of a fractional derivative of order a > (e.g. [22]) that it yields zero when 
applied to a constant. This property requires that 



/. 



" f{x)dx = Const X a-"" (106) 



which is obviously not solved for a > by any real valued power-function f{x) = const x~°'~^ G R since its 
primitive in (106) diverges at x = 0. In contrast ba{x) solves (106) 



RelTl ^^r^l=^4^7 WlS°l=0, a>0 (107) 

for any e > since the lower limit of the integral is purely imaginary and the upper limit vanishing for a > 0. 
Hence (95) defines for a > an integration kernel for a fractional derivative which yields zero when applied to 
a constant which is due to the vanishing of (99) which we can indeed verify directly. We analyze in appendix 
8.1 further properties of ba- 

Let us now return to the functions qn{x) of (86): Taking into account (95) have with a = nS > —1 and 



with 



bnsix) = -r^Six) = lim Re I — } , nS > -I (108) 

where Ag is given in (24). For n = 1 we obtain 

1 5' f 1 

bsix) = -rC,5{x) = - lim Re I -—. —itt \ (109) 

being in accordance with qi{x) = C5{x) = —/S.(^s,h,c,)^{x) obtained by evaluation of (7) or (8) for x 7^ 0. Prom 
(108) one can recover the static Green's function by (13) 

g-i(x) = = L-H{x) = ^i^e T, , \, , \ , < ,5 < 2 (110) 

^ ' As ^ ' -kAs \(|x| +ze)i-'5j V ) 

which indeed recovers expression (27) for this Green's function q-\{x) = g{x). (110) takes into account in a 
"correct" way the singular behavior at x = and is therefore more general as (27) holding only for x 7^ 0. 
Whereas the 6^5 (x) hold for nS > —1, the functions qn{x) = A'^bns{x) are only physically meaningful where 
As > and where (7) exists namely < S < 2. The range of validity of the 6q,(x) can be mathematically 
further extended to all a < — 1 if one utilizes instead of (11) the regularized definition of the F-function (34). 

The performed analysis gives reason to the following observation: There is a general link between non- 
integer a including odd integers a and the non-locality and self-similarity of physical properties. This behavior 
is also refiected by the non-locality of our Laplacian (7). Intuitively in a wider sense one is tempted to state: 
even integer-exponents a (case (iii)) correspond to local ("regular"), non-integer and odd exponents a (case 
(iv)) to non-local, self-similar ( "irregular" )^^ physical behavior. 



The meaning of the term "irregular" is here used in the sense of "erratic", it does not mean the absence of rules. 
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7 Conclusions 



We have analyzed an elastic medium with self-similar elastic energy-density. The source of the physics in 
this continuum is the self-similar Laplacian (7) having the form of a non-local convolution with a power-law 
convolution kernel. This Laplacian can be represented in terms of fractional derivatives and is obtained by the 
continuous limiting case of the discrete self-similar Laplacian introduced in [13]. The resulting equations of 
motion are sclf-adjoint partial fractional differential equations which are only in the range < 5 < 2 physically 
meaningful. 

We deduced closed form solutions (i) for the static Green's function (displacement field due to a unit 
(5- force), (ii) for the Cauchy problem (two convolution kernels "propagators"), (iii) for the solution of a self- 
similar diffusion problem. The solution was obtained in terms of a convolution kernel solving the initial value 
problem. It has been found that the obtained convolution kernel represents normalized Levi-stable probablity 
distributions with inifinite variances (mean fluctuations) in their range of deflnition < 6 < 2. 

As a spin-ofF result we obtained a set of functions bnsix) which we denoted as self-similar potentials 
appearing as regularized Fourier transforms of \k\°'. These functions arc to be conceived as generalized 
functions in the distributional sense of Gel'fand and Shilov [23]. Self-similar potentials are special solutions 
of the self-similar Poisson's equation and play in statics the analogous role as Newtonian potentials for the 
traditional Poisson equation. 

It has been demonstrated that the concept of self-similar functions introduced in [13] is useful to describe 
material systems with self-similar interactions. Our approach opens the door to a continuum description of 
a rich fond of new physical especially dynamic phenomena in such systems. Especially interesting arc the 
propability-distributions obtained as solutions of the self-similar (fractional) diffusion equation (68). These 
probability-distributions are Levi-stable which are known to be a characteristic feature of the statistics of 
irregular and erratic motions characterized by discontinuous trajectories (Lcvi-fiights) [1, 3]. It is mention- 
worthy that our diffussion model recovers the range < 5 < 2 of admissible Levi parameter 6 given previously 
in the literature (e.g. [1, 3]). 

We have further formulated the ID continuity equation for the self-similar diffusion problem and formu- 
lated the self-similar probability-flux (Eq. (71)). An interesting open point which can be treated by the 
present approach remains the formulation of the non-local self-similar Gauss-theorem in dimensions higher 
than one, and more general to formulate continuity equations and related problems in self-similar continua 
embedded into the 2D or 3D space. 

The present approach can be a point of departure to answer such questions in a rather rigorous way. 
For example the present approach can be considered as a point of departure to any general self-similar field 
theory such as continuum mechanics and electrodynamics. We hope the present paper stimulates work in 
such directions. 



8 Appendix 

8.1 Fractional derivative 

In view of the analysis of the last section we can easily recall a fractional derivative D^f{x) which fulfills for 
a > all desirable properties of a derivative. We call positive orders (fractional) derivatives and negative 
orders (fractional) integrals. Let us consider the case a > — 1 which can however be extended to a < —1 in 
the spirit of generalized functions [23] as outlined in appendix 8.4. 
Let us first of all determine 

ya{x) = D^S{x) , a>-l (111) 

where we call yaix) a fractional derivative. For a = the solution of should be ya=oix) = S{x). 

From (111) follows then for the deflnition of this fractional derivative applied on a function f{x) the 
convolution 

/oo 
ya{x-T)f{T)dT, a>-l (112) 

-oo 

with the kernel ya{x) (115) to be determined. Further we assume that D^e^^^ = {ik)°'e^^^ so we can write 
yaix) write as a Fourier integral 
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which we can rewrite 



ya{x) = — / e'^^itkrdk (113) 
y^{x) = -1 hm r e-'^k"^ ( e^(*^^+'^) + e-'^^'^+^A dk (114) 

2tt e^0+ Jo ^ ' 



which can be further written as 



Vaix) = — hm Re \ . \ , a > -1 (115) 



For X y^O (115) assumes the form 



ya{x) = , sinTra (116) 

7^ a+i 



We verify that for integer a = 0, 1, 2, .. G Wo (116) is vanishing for x 7^ and (115) yields 

yjx) = - — hm Re < -—^ )■ = hm Re < r )■ = - — 6(x) (117) 

reproducing in the integer-case the correct locahzed the integer-order derivatives. 
8.2 Convergence of series (50), (51) 

It stih remains to verify the convergence of the series (50) and (51): To this end we consider series (51) which 
has the form 



|ofe*)=*w-^^™ (118) 



Im \ P I 

7r\x\ I \ 

where ^0(2^) = S(x) and P is only a function of ^ = ^^jl^) where P is given by 

\x\ 



^(0 = f:|^r, 0<S<2 (119) 



In order to judge the convergence of this series the asymptotic representation of for n >> 1 wih be 



(nsy. 

useful. To this end we put s{n) = 2n(l — |)) with nS = 2n — s 

{nSy. (2n-s(n))! 1 



(2n)! 2n! (2n)2"(i-f) 



6\ 5 



< (5 < 2 (120) 



Prom (120) we can determine the radius of convergence pc of the series (119) where the series OnC" 
converges absolutely for |^| < pc with 

Pc = hm 1-^1 = lim (2ne)^"5 = 00 , 0<6 <2 (121) 

since 1 — | > 0. Hence series (119) converges for all ^ absolutely which guarantees the convergence of (50) 
and (51) for all x and t for 6 being in the interval < 5 < 2. 
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8.3 Convergence of series (85) 

We consider the convergence of (85) 



Wix,t) = J2i-ir-qnix), 0<d<l (122) 

n=0 



which we can rewrite in the form 



TTld 

Wix, t) = Six) - -L-Im{w{^^)} (123) 



7vi6 

where w{^) with ^ = "^^^^^ has the representation 



^(0 



^ (ndV 

^(-l)"^r, 0<5<1 (124) 



n=l 



The radius of convergence rc of this series is 

lim M! + 1)' = + 1) ^ oo , 0<S<1 (125) 

n-s-oo n! ((n + l)(5)! ^ ^((n + l)5)! ^ ^ 

Hence (124) converges for all ^ absolutely uniformly. Hence (85) converges in the entire domain — oo < 
X < oo and t > for 5 being in the interval < S < 1. 

8.4 Some remarks on the Fourier-integral (94) for a < —1 

A further observation is important in ba{x) of (94) which is the Fourier transform of exists for a > 
— 1. For exponents a < —1 the integral (94) diverges. However, we can define a Fourier transform of 
Re{{—i)"{\k\ + ie)"} in the limiting case e — > 0+ which exists also for a < — 1. 

We generalize the definition of ba{x) (94) for a < — 1 by replacing — )• ^^ea Re{e — i\k\)'^) in the form 

baix)= lim 3_— i?e((-i)" / (k + ie)" cos kx dk] , a G IR ^ ±1, ±3, .. ± (2n + 1), .. (126) 

e->0+ TT COS -g- K Jo J 

which exists a G IR where the zeros of the cosine have to be excluded. For a > —1 (126) is identical with (94) 
and for a < —1 it complements the representation (94) which diverges in this case. Let us especially evaluate 
(126) in the "forbidden" case a < —1: Taking into account 

(e - ilkl)'^ = — / e-^(^-^l*^l)r-°-MT (127) 

r(-Q;) Jo 

existing for a < —1 we obtain the representation 

bJx)= lim — rRe\ dre-'V-^-M cos A;x e'^^dA; ^ (128) 

^ e_,o+7rcos^r(-a) [Jo Jo J ^ ^ 

which we evaluate by using 

Re / cos kx e'''^dk = -5{t - \x\) (129) 
Jo 2 

and arrive at 

^«(^) = 2col^ r(-a)' ' a<-l,^ -1, -3, .. - (2n + 1), .. (130) 

where r(— a) is for a < —1 well defined in (11). 
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By taking into account the Euler relation (34) we arrive at 



ba{x) = - 



a! 




a<-l,^-l,-3,.. -(2n + l),.. 



(131) 



a+l 



which has the same form as expression (98) which holds for a > —1 but with the generalized definition of a\ 
given in (34). In conclusion (131) holds for a G IR ^ ibl,±3, .. ± (2n + 1), ... For cases a being even integers, 
expression (131) is not valid, we refer in this context further to [23] where the Fourier transforms l/|fep"'"'"^ 
(n = 0, 1, 2..) are further analyzed. 
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